Correlation length by measuring empty space in simulated aggregates 
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We examine the geometry of the spaces between particles in diffusion-limited cluster aggregation, 
a numerical model of aggregating suspensions. Computing the distribution of distances from each 
point to the nearest particle, we show that it has a scaled form independent of the concentration 
<f>, for both two- (2D) and three-dimensional (3D) model gels at low <f>. The mean remoteness is 
proportional to the density-density correlation length of the gel, £ , allowing a more precise measure- 
ment of £ than by other methods. A simple analytical form for the scaled remoteness distribution is 
developed, highlighting the geometrical information content of the data. We show that the second 
moment of the distribution gives a useful estimate of the permeability of porous media. 

PACS: 61.43.Hv, 82.70.Gg, 47.55.Mh, 02.70.-c 



Substantial attention has been paid to cluster-cluster 
aggregation as a computational model of gelation in 
colloidal suspensions, aerogels, aggregating emulsions, 
etc. pi. There is a large literature analysing the struc- 
ture of individual aggregates, in terms of fractal geometry 
^,|). However it is now properly understood that, 
in irreversibly aggregating systems at non-zero particle 
concentration 4>, clusters, while fractal at intermediate 
length scales, must inevitably become space-filling at the 
largest scales, as clusters pack to form a 'gel' with a well 
defined correlation length £. Direct analysis of the gel 
structure remains rare |5(], although it is fundamental to 
a gel's macroscopic physical properties. Colloidal gels 
and aerogels are porous media and, as such, the geome- 
try of their 'pores' or vacancies is of prime importance, 
influencing the flow of permeating solvents and gel sta- 
bility under external stresses such as gravity raffi]. Here, 
using the diffusion-limited cluster aggregation (DLCA) 
model we examine the structure of the space inside 
the porous gel. DLCA is perhaps the simplest realistic 
model of aggregation in colloidal suspensions, and has 
been shown to qualitativelyreproduce many features of 
experimental systems We show that examina- 

tion of the void structure in the simulated gel leads to 
a measure of the correlation length £ of unprecedented 
precision, while the scaling with <j> also gives an estimate 
of the cluster fractal dimension at scales smaller than £. 
We further find that, at low (/>, the pore structure is in- 
dependent of <fi, i.e. DLCA particle gels have statistically 
identical structure regardless of particle concentration. 

In fractal aggregates, there is no single well-defined 
pore size, since spaces of all sizes are present. To inves- 
tigate the distribution of void sizes we require a suitable 
measure; but it is difficult to avoid some arbitrariness 
in delineating the extent of a void pj] , p^ |. Marine car- 
tographers face a similar problem: where does one sea 
end and the next begin? We bypass the problem as fol- 
lows. If an adventurer adrift in a boat sights land nearby, 
(s)he cannot be sure whether (s)he is on a small lake, or 
near the edge of a large ocean. However, if (s)he knows 
that the nearest land is far away, (s)he can infer the ex- 



istence of a large body of water. Hence information on 
the distribution of void sizes can be gained by measuring 
the remoteness, i.e. the distance to the nearest particle, 
of all points in an aggregated system. This distribution 
of remoteness should not be interpreted, as in Ref . [p"3| , 
as the distribution of void sizes. The former is only an 
indirect measure of the latter. 

The DLCA algorithm is implemented as follows. A 
square (2D) or cubic (3D) grid, with periodic bound- 
ary conditions, is occupied by a randomized population 
of particles at concentration <j>. Each particle performs 
a random walk on the lattice until it finds a particle 
on a neighbouring site. The two then stick irreversibly. 
A cluster of particles thus formed continues its random 
walk, at a rate inversely proportional to its radius of gy- 
ration, forming ever larger aggregates on contact with 
neighbouring clusters. The growing floes have a frac- 
tal geometry characterised by a dimension df which is 
lower than the space dimension d, so that, on the scale 
of a cluster of 'mass' M and radius r, the mean con- 
centration (~ M/r d ~ r d, ~ d ) decreases as the cluster 
grows. When the clusters grow to such a size (£) that 
their mean concentration (~ £ d f~ d ) matches the overall 
system concentration (f>, mass conservation demands that 
the clusters must fill space. Assuming a concentration- 
independent fractal dimension (evidenced in the limit of 
low 4> ||[Hj|]), the correlation length in the 'gels' varies 
as £ oc <fr^i To avoid artifacts in the simulations, 

the correlation length £ must be both significantly larger 
than the lattice parameter a, and smaller p3| than the 
system size L, L ^> £ ^> a. 

We analyse the 'gel' structure after all particles have 
aggregated. The fact that the particle positions are con- 
fined to the lattice is irrelevant when finding the distribu- 
tion of remoteness; we visit all points in the system, and 
measure the Euclidean distance to the nearest particle 
(contrast Ref. [Of). For each gel, we compile a histogram 
of remoteness MR), which is normalized so that h(R) dR 
is the probability that a point chosen at random has a 
remoteness in the range R ± dR/2. It is significant that 
we find h(R) to be normalizable. This is not the case for 
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similar analyses of isolated fractal clusters [ |ll| , |i^ , |l4j , 
where an artificial boundary must be imposed as the 
empty space surrounding a cluster is, in principle, lim- 
itless. The distributions we measure contain a natural 
cut-off, revealing that no points in the system are more 
remote than some finite distance. 

Given the hypothesis that the correlation length is the 
only relevant length scale present, the cut-off distance 
must be a fixed fraction of £, as must other lengths such 
as the mean remoteness (R). To check whether a single 
length scale characterises the geometry at various (f>, we 
re-scale each histogram by its mean (R), and define a 
scaled remoteness distribution p(x) with unit mean and 
unit norm, by p(x) = (R) h((R)x), where x is remoteness 
measured in units of the mean (Figs. and || for d = 2 
and d = 3 respectively). Points at the smallest values of 
x are affected by grid artifacts, but otherwise, the data 
collapse very well onto a single master curve for each 
value of d. 




1 2 3 4 5 

FIG. 1. For simulations in d = 2, the distributions, p(x), 
of remoteness in units of the mean remoteness, for concentra- 
tions 4> = 0.10 □, 0.15 O, 0.17 A, 0.20 <, 0.23 v, 0.25 >, 0.27 +, 
0.30 O- Solid line: one-parameter theoretical fit. Inset: The 
same graph on a log scale and over a wider range, making the 
tail of the distribution more visible. 

No fitting was required to achieve the data collapse 
m Figs. and |, since the scale factor (R) is a well de- 
fined property of the gel. Furthermore, when we plot this 
length scale, the mean remoteness, against concentration 
cj) in Fig. [|, we obtain remarkably clear power-law be- 
haviour. As discussed, a power-law relationship between 
4> and any length scale proportional to £ is expected at 
these low concentrations, but the almost complete ab- 
sence of noise was quite unforeseen. Compare equiva- 
lent plots in the literature || using the length scale ob- 
tained from the peak in the structure factor. We deduce 
that moments of the remoteness distribution are precise, 
noise-free measures of the length scale present in the gel. 



Further confirmation that (R) is a measure of, though 
much smaller than, the correlation length £ comes from 
the gradients of the lines in Fig. [|| which are respectively 
-1.67 ± 0.05 and -0.73 ± 0.03 for d = 2 and d = 3. 
Given our expected gradient —l/(d — df), this implies a 
fractal dimension df — 1.40±0.02 for d = 2, in agreement 
with the accepted value of 1.42 (see e.g. |Q), and df — 
1.63 ± 0.05 for d = 3, marginally below the accepted 
zero-concentration value of 1.80. 
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To find the constant of proportionality relating (R) 
to the more commonly used measure of the correlation 
length, from the peak in the structure factor at wavevec- 
tor <? max , we compute the ratio (R)/£ = {R)q max /2ir for 
each gel, and find a constant value 0.034±0.004 for d = 2, 
and 0.121 ±0.007 for d = 3. We also measure higher mo- 
ments of the remoteness distribution (table |) . 

To understand the shapes of the master curves in 
Figs. and 0, we take as our starting point the fact that 
the gelled aggregate is self-similar on scales smaller than 
the correlation length. On some scale, we may picture the 
structure as a set of amorphous blobs and cavities. These 
cavities contain points with a range of remoteness, char- 
acterised by a normalized remoteness distribution f(R). 
The precise form of f(R) depends on the shape of the 
cavity. If we l zoom ivH to a smaller length scale <r, we no- 
tice that the blobs are themselves aggregates of blobs and 
cavities. By self-similarity, those smaller cavities have a 
normalized remoteness distribution f(R/a)/a, i.e. of the 
same shape as the larger cavities, but scaled down. Thus, 
compiling the overall histogram of remoteness, h(R), we 
get contributions proportional to f(R/a) from all scales 
<7, up to a largest scale x ~ £■ The magnitude of the con- 
tribution from each length scale is as yet undetermined. 
Let us call it c(<j). Hence, 
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h(R) 



c(a) f(R/a) da. 



(1) 



We have expressed h(R), the remoteness distribution for 
a complex fractal aggregate, in terms of f{R), the dis- 
tribution of remoteness in some 'primal' cavity with a 
simple (albeit unknown), non- fractal geometry. Equa- 
tion (0) may be regarded as a definition of f(R). 

Since the geometry is self-similar on all length scales 
a in the integration, < a < x (ignoring a lower cut-off 
due to the grid), c(er) must be a function with no charac- 
teristic length scale, i.e. a power of a up to the cut-off x- 
In fact, we have no prior justification for the cut-off to 
be sharp since, in principle, scale-invariance could be lost 
gradually over a range of scales up to x- Nevertheless, 
the hypothesis will be validated by the data. We define 
the as yet undetermined exponent v by c(a) = cqo~ v ~ 2 . 
The constant coefficient cq is fixed by the normalizations 
of h(R) and f(R) as Co = vx v - As mentioned, the 
distribution of remoteness h(R) is not synonymous with 
the distribution of void sizes. We may identify c(cr) as 
the latter, while f(R) contains information on the void 
shapes. 

Recasting (nf) in terms of the scaled master curve 
p(x) = (R) h((R)x), and putting z = (R)x/a, yields 



p( x ) = - ( — 
x \x 



/xo 



(2) 



where x$ = x/{R), which is not a free parameter: from 
Eq. (|l|), Xq 1 = {x)fv/{v + 1) where (x) j is the mean 
of the distribution f(x). Without loss of generality, we 
set the size scale of the cavity described by f(x) so that 




0.02 0.1 
FIG. 3. Log-log plot of {R)/a vs. (f> for 2D (upper set) and 
3D (lower set) simulations. 



d 


(R 2 ) 1/2 /(R) 


{R A ) 1/A / (R) 


(R 4 ) 1/4 /(R) 


(R 5 ) 1/5 / (R) 


2 


1.42(3) 


1.77(5) 


2.07(8) 


2.35(10) 


3 


1.141(9) 


1.25(2) 


1.35(3) 


1.44(3) 



TABLE I. Ratios of the moments of the remoteness distri- 
bution, averaged over all simulations (error in last digit). 



the most remote point within it is separated by unity 
from the nearest boundary. Hence f(x) — for x > 1, 
implying p(x) — for x > xn. So we interpret xq as the 
point on the cc-axis of Figs. and | at which the graph 
falls to zero, i.e. the value, in units of (i?), of the natural 
upper cut-off in the remoteness distribution. 

The value of the parameter v can be calculated by 
examining the small- a; behaviour of Eq. (Q). An ob- 
ject closely related to the remoteness distribution is the 
Minkowsky cover [|l2],[l4 15 1, the locus of points lying 
within a fixed distance, R say, of a fractal set. For 
small R, the volume of the Minkowsky cover grows as 
Rd-df ||Ji4||rj]. Our distribution h{R) is proportional 
to the perimeter of the Minkowsky cover, i.e. the deriva- 
tive of its volume with respect to R. Hence, we require 
p(x) oc x d ~ d f~ 1 in the limit of small x. As x — > 0, the 
integral in Eq. (^), assuming it is convergent, tends to a 
constant, so that p(x) oc x 1 '^ 1 and we identify v = d — df. 

It remains only to determine the form of the primal 
distribution f(x), i.e. the geometry of the basic cavity. 
We know it is normalized and must vanish for x > 1 . We 
shall discuss some further features but, because Eq. (Q) 
already embodies many of the correct properties for a 
gelled aggregate, we find the results to be insensitive to 
the precise form of the function /. 

Let us examine how f(x) vanishes as x approaches 
unity from below. Consider a circular cavity with unit 
radius in 2D, inside which are drawn contours of equal 
distance from the edge, i.e. concentric circles. The cav- 
ity's remoteness distribution f(R) is proportional to the 
length of the contour at a distance R from the edge. 
Hence, as we approach the centre of the cavity, R — > 1, 
f(R) falls linearly to zero. This would also be the case, 
in 2D, for a triangular cavity, or any regular polygon. 
In the centre (most remote part) of a 3D cavity, on the 
other hand, the distribution vanishes quadratically, since 
it is proportional to the area of an equal-remoteness sur- 
face. Hence, one might consider substituting the nor- 
malized distribution /trial = (1 — x) d ~ 1 d for < x < 1 
into Eq. (||), thus building the fractal texture from a set 
of regular 'lagoons' of all sizes. This formula was not 
used for the curves in Figs. |l| and |[ despite good agree- 
ment with the 2D data, as it has shortcomings in 3D. 
For d = 3, v » 1.2 so that, as x — > 0, the integral in 
Eq. (||) diverges if /(0) is finite. At small x, we demand 
p(x) ~ x d ~ df ~ 1 , which requires p to be an increasing 
function of x for d = 3, where df « 1.8. When we con- 
struct a surface of constant distance around a tree-like 
'branch', near the branch the area of the surface must 
grow as the distance from the branch is increased, be- 
cause the branch is convex. But the negative gradient of 
/trial describes only concave lagoons. According to our 
data, such a function builds an adequate picture of the 
spaces in a 2D aggregate, since the aggregate with df > 1 
is sufficient to partition 2D space into concave lagoons of 
all sizes. But in 3D, branches form an imperfect cage 
which has significant convex structure. 

To acknowledge the existence of both concave and con- 
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vex parts of the aggregate, we use the simplest possible 
function f(x) that vanishes in the correct way (with a 
power d — 1) both at x = and x = 1. We construct a 
piecewise smooth function that is the lower envelope of 
two polynomials: c\x d ~ x and C2(l — x)^ 1 with constant 
coefficients c\ — j3 1 ~ d d and C2 = (1 — f3) 1 ^ d d chosen to 
ensure that it is continuous and normalized. The posi- 
tion of the cusp where the two parts meet, at x = j3, is 
the only free parameter of our model. The mean of this 
function is at (x) f — (l+[d— l]0)/(d+l) so that we may 
eliminate the parameter (3 in favour of Xq which we fit to 
the measured value of the cutoff in the master curves. 

We find best values xq = 6.5 ± 1.0 and 3.5 ± 0.2 for 
d = 2 and 3 respectively. These one-parameter fits are 
the curves shown in Figs. and |. Our single formula 
reproduces the data very accurately both for d = 2 and 
d = 3, even in the tails of the distributions, until the data 
become noisy. 

The second moment of h(R) has a special significance 
to porous media. At low Reynolds' number, the mean 
flow rate, (it), of a solvent of viscosity r/ through a porous 
medium subject to a pressure drop per unit length VP, is 
given by Darcy's law |16|, (u) = kVP/rj, where the per- 
meability k is a squared length scale that depends only on 
the geometry of the porous medium. To avoid complex 
hydrodynamic calculations, k is often estimated by the 
Carman-Kozeny (CK) relation @, k » (1 - </>) 3 /5S 2 , 
in terms of the specific surface area S. However, since 
the flow velocity at any given point is not as dependent 
on the total surface present as it is on the distance to 
the nearest obstacle, the CK relation mis-calculates the 
porosity in materials with more than one characteristic 
pore size. For instance, the flow rate through a set of 
parallel cylindrical holes, of various radii a, can be cal- 
culated exactly using Poiseuille's formula yielding 
k = (1 — 0)(a 4 )/8(a 2 }, which, expressed in terms of the 
remoteness, is k = 3(i? 2 )/4. By substituting the sur- 
face area of the holes into the CK formula, one obtains 
kcK ~ (1 — </>)(a 2 }/10, which approximates the exact 
value only for a sharply peaked distribution of radii, for 
which (a 4 ) « (a 2 ) 2 - The approximation fails too for a 
polydisperse set of parallel-sided channels, where the ex- 
act result can be expressed as k = (R 2 )/2. We there- 
fore propose an approximate rule of thumb k ss 0.6(i? 2 }, 
which continues to hold for polydisperse pores and for 
rough surfaces while, in both cases, the CK equation fails. 
The new formula has been shown to yield results of the 
correct order of magnitude for fractal particulate gels || . 
However, like the CK equation, the new formula lacks in- 
formation on pore topology, so must also fail in certain 
cases, e.g. where some pores are closed. There is then no 
alternative to a full hydrodynamic computation. 

The shape of the remoteness distribution carries infor- 
mation on the geometry in the DLCA simulations, and 
confirms (by scaling onto a master curve) that, statisti- 
cally speaking, all the simulations at the low concentra- 
tions considered, give rise to the same geometry: the par- 
ticle gel 'pore' structure is independent of concentration. 



The mean remoteness, an easy quantity to measure in 
the simulations, is proportional to the correlation length 
(at least in the range of concentrations studied), hence 
giving a measure of the gel's characteristic length scale 
that is considerably less noisy, and perhaps more easily 
determined, than other measures. The second moment of 
the remoteness distribution provides an estimate of the 
permeability of the gel. These coordinate-independent 
measures of real-space length scales could be applied to 
a range of other phenomena, such as the void geometry 
in reaction-Yamied aggregation |19| and in earlier stages 
of aggregation ||], or the ^-dependence at higher con- 
centrations H when the structure is no longer controlled 
solely by the zero-concentration critical phenomena. 
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and the Royal Society of Edinburgh. 
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